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Abstract—Quality control is critical for all industrial processes, 
but often times defect detection is labor intensive. A novel 
approach to industrial defect detection is proposed using a Digital 
Noise Radar (DNR), coupled with Radio Frequency Distinct 
Native Attribute (RF-DNA) fingerprinting processing algorithms 
to non-destructively interrogate microwave devices. The DNR is 
uniquely suitable since it uses an Ultra Wideband noise waveform 
as an active interrogation method that will not cause damage 
to sensitive microwave components. Additionally, it has been 
demonstrated that multiple DNRs can operate simultaneously in 
close proximity, allowing for significant parallelization of defect 
detection systems resulting in increased process throughput. 
Using this method, 100% sampling for quality control may be 
attainable in many cases. 

The ability to classify defective units from properly functioning 
units was demonstrated in [1] with potential applications includ- 
ing assembly line inspection of automotive collision avoidance 
systems, wireless or cellular antenna defect detection during 
manufacture, and phased array element defect detection prior 
to RF system assembly. However, prior research into active 
interrogation has been strictly empirical. This paper will develop 
an analytical model and simulation for interrogating a passively 
terminated antenna with an overall objective of improving clas- 
sification performance through optimization of the interrogation 
signal bandwidth. 

Keywords: Noise Radar, Device Classification, RF-DNA, Ultra- 
Wideband 


I. INTRODUCTION 


RF-DNA fingerprint processing is a method to extract 
features received from a Radio Frequency (RF) signal [2]. 
The concept is based on the fact that every device emits 
signals that have unique characteristics that can be used to 
distinguish that device from other, similar devices [3]. The 
term “fingerprinting” is used because the concept is similar 
to human fingerprinting, where a person’s fingerprints can be 
used to distinguish them from another person [2]. RF-DNA 
fingerprinting has been in development at Air Force Institute 
of Technology (AFIT) since 2006 successfully demonstrating 
feature extraction from RF emissions received from various de- 
vices including IEEE 802.11 WiFi, IEEE 80211.15 BlueTooth, 
IEEE 802.16 WiMAX, Global Systems for Mobile Communi- 
cations (GSM) cell phones, and Radio Frequency Identification 
(RFID) emitters [4]. To date, all of these methods have been 
performed in a passive manner, whereby the signals analyzed 
by the RF-DNA fingerprint processing algorithms were already 
generated by the device under normal operation. This works 








well for devices that normally broadcast some RF emission 
that can be exploited. However, a whole class of devices 
exist that do not normally emit RF radiation, such as properly 
shielded RF amplifiers, mixers, oscillators and filters that make 
up a typical receiver front-end. The potential exists for these 
devices to be classified using an “active” interrogation method, 
such as a radar transmission. The RF energy of the radar will 
be reflected by the device in a manner that can be exploited 
by RF-DNA fingerprint analysis [1]. 

Instead of a traditional narrow-band radar signal for inter- 
rogation, this paper investigates the use of an Ultra Wideband 
(UWB) noise signal from a DNR as the interrogation wave- 
form, specifically the effect of bandwidth on classification 
performance. An UWB signal has a much larger frequency 
content than a traditional narrow band signal. As such an UWB 
signal may provide increased classification performance as an 
interrogation signal because of the larger instantaneous band- 
width [3]. For instance, the UWB signal may excite resonant 
structures in a specific target (such as a Log Periodic Antenna 
(LPA)), resulting in larger signal energy returns within specific 
frequency bands which can be used for classification of the 
radar target [5]. The bandwidth requirements of an UWB 
Digital Noise Radar drive the design requirements of the 
system. Ideally, a noise source with infinite bandwidth with 
a flat spectral response is available for classification, but this 
is unrealizable; all signals are band limited and it is important 
to be able to thoroughly characterize the impact of signal 
bandwidth on classification performance. 


Il. RF-DNA FINGERPRINTING 


Statistical RF-DNA fingerprints consist of features that are 
generated by the statistical properties of fixed regions within 
the received signal [2]. The RF-DNA fingerprint process 
extracts features using three instantaneous responses: the In- 
stantaneous Amplitude (IA) response, a(n), the Instantaneous 
Phase (IP) response, ¢(n), and the Instantaneous Frequency 
(IF) response, f(n). The standard deviation (ø), variance (o?), 
skewness (y) and kurtosis (x) are generated for each response. 
These statistics are calculated using the following equations: 





n)—p), (1) 
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where e(n) is the normalized data sequence with N, samples, 
and each of the above statistics is calculated over an equal, 
contiguous, sub-region of e. These sub-regions are deter- 
mined empirically and are application specific. For the UWB 
experiment in section V five equal sub-regions were used. In 
each of the five sub-regions, standard deviation and the three 
statistics calculated in (1), (2) and (3) are concatenated to form 
“distinct native attribute marker” to the sub-regional RF-DNA 
vector, i, according to [2]: 
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After which, each distinctive native attribute marker vector 
is concatenated to form a “composite characteristic vector” 
according to [2]: 


Fe = |F Fr, i++ Fry | : 
Ayia Ryat 1z4(Nr+1) 


in which C denotes the analyzed response (IA, IP or IF), 
and Np is the total number of sub-regions. Each of the 
composite characteristic vectors are finally concatenated into 
one statistical fingerprint (of 60 elements for 5 sub-regions) 
according to [2]: 


F= [r° Be: Fi, (4) 


A. Multiple Discriminant Analysis / Maximum Likelihood 


Multiple Discriminant Analysis (MDA) is not, by itself, 
a classifier [6]. It is an algorithm to achieve dimension- 
ality reduction of a multi-dimensional dataset such that a 
more traditional classifier can be applied, e.g., a Maximum 
Likelihood (ML) classifier. Multiple Discriminant Analysis / 
Maximum Likelihood (MDA/ML) aims to project from a d- 
dimensional space to a (c — 1) dimensional space, where c 
is the number of distinct classes (it is assumed the d > c, 
or some space of dimension greater than c). Figure 1 is a 
example of MDA/ML applied to a 3-dimensional dataset. In 
the figure, the MDA/ML algorithm is used to project the 
3-dimensional space into two, 2-dimensional sub-spaces. The 
2-dimensional sub-spaces are defined by their norm vectors, 
W, and W3 respectively. The algorithm attempts to find the 
optimal sub-space for classification, by finding one with the 
greatest separation of the projected distributions of the classes. 
For instance, discrimination along the W, sub-space would be 
considerably easier than across the W2 sub-space that still has 
overlapping projections. For higher class problems, it becomes 
difficult to visualize in a 3-dimensional space. 





Fig. 1. MDA/ML classification of fingerprint in a 3-class problem [6] 
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Fig. 2. Proposed transmission line model for target antenna under test 


II. PROPOSED TRANSMISSION LINE MODEL 


The “target” for the research conducted is a passively ter- 
minated UWB antenna. Previous experimentation using both 
LPA and Ultra High Frequency (UHF) antennas as the transmit 
and receive antennas on a UWB Digital Noise Radar has been 
conducted at AFIT for several years [3], [5], [7]-[10]. For 
this model, an electromagnetic plane wave signal, generated 
from an UWB digital noise radar, propagating in free-space, is 
incident on the bore sight direction of a passively terminated 
antenna. The reflected signal received at the noise radar is 
then processed using a digital correlator which can be used for 
various purposes to include ranging information and RF-DNA 
based classification. 


A. Scatter 


The use of an antennas as a target model effectively cat- 
egorizes the system response as an electromagnetic scatter 
problem. Scattering is the effect of removing energy from 
a propagating electromagnetic wave, and re-radiating it with 
potential changes in direction, phase or wavelength [11]. For 
electromagnetic scattering, some of the energy incident on 
the antenna is absorbed and some is re-radiated towards the 
receiving antenna. The three major types of scatter from the 
antenna are: antenna-mode scattering, as, antenna-structural 
scattering, ca and RCS grating [11]. The structural mode 
scattering is from the physical shape of the antenna, material 
composition, geometry as well as the orientation. The antenna- 
mode scattering is due to the internal reflections inside the 
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antenna, including reflections off of the termination for the 
proposed passively terminated antenna model. Grating lobes 
are areas of increased radiation in non-bore sight directions 
on linear arrays, and have no real relevance for the proposed 
system model. Additionally, it is assumed that the incident 
wave polarization is matched to the “target” antenna, since 
some antennas (such as dipoles) have little to no backscatter 
for unmatched polarizations [11]. 

For the structural scattering mode, which is based on 
various structural characteristics of the antenna, the Radar 
Cross Section (RCS) will remain fixed at a given incident 
angle. However, for the antenna mode scattering, the RCS is 
based on the termination load and the antenna impedance. For 
the purposes of this research, it is assumed that the antenna 
impedance is designed to be 50Q across the bandwidth of the 
UWB antenna. 


B. Filter Model 


The spatial aspect of the filter model concept concerns 
antenna directivity. For the research conducted within this 
paper, the incident electrical field is modeled as a plane wave 
incident on the bore sight of the antenna, or the axis with 
maximum gain. A key assumption is that the direction of 
the bore sight does not change with frequency. Therefore, 
both the transmit mode and receive mode of the antenna can 
be modeled with an impulse response h(t, 0, ọ%) in the time 
domain and corresponding transfer function H (w, 0, %) in the 
frequency domain. The impulse response and transfer function 
is reciprocal for both transmit and receive (Hrz(w,0, Y) = 
Arz(w,9,~)), with the exception of the signal vector direc- 
tion. Since we are only concerned with the boresight direction 
of the same antenna, separate transmit and receive transfer 
functions are unnecessary [12]. 

The transfer function of the antenna only models the the 
energy propagating through the antenna. The received transfer 
function, Hp, is defined as the relationship between the inci- 
dent electromagnetic field, Erz(f,r) to the voltage output at 
the antenna terminals, U pz (f), and likewise for the transmitted 
transfer function, Hy, [12]. These transfer functions do not 
take into account the structural-mode scatter or the antenna- 
mode scatter mentioned previously. However, for the proposed 
antenna-target, antenna-mode reflections off of the antenna 
loading will be captured by the cascaded receive and transmit 
transfer functions. 

In order to investigate the effect of signal bandwidth on 
classification capability of the RF-DNA tools, the antenna can 
either be modeled as a low-pass filter or a band-pass filter. 
[5] showed that the bandwidth of the LPA antennas used in 
the AFIT DNR, and proposed as the target antenna, has an 
approximate bandwidth of 400 - 750 MHz, which can be 
modeled as a bandpass filter. The cross-correlation response 
is affected by the center frequency of the interrogation signal, 
therefore the bandpass filter can be tuned not only to a different 
bandwidth, but also a different center frequency. Finally, to 
simulate the effect of a more ideal UWB signal with a much 
larger frequency content, a low pass filter can be employed. 








The resonant modes of a LPA will be difficult to model with 
a bandpass filter, but for the purposes of this research, the 
resonant effect is not really necessary and could serve only to 
complicate the model. 


C. Time Dispersion 


Modeling an antenna as a filter means that it has a mag- 
nitude and phase response, H(f) = A(w)e/*) [13]. Filters 
will distort the input signal, instilling both amplitude and phase 
distortions. These distortions comprise the group delay of the 
filter, which indicates how much time delay the filter imposes 
on the input signal at specific frequencies. This is a good 
measure of the time dispersion of the antenna. A linear filter 
will have an ideally constant group delay, meaning that all fre- 
quencies are time delayed by approximately the same amount 
resulting in little distortion of the output signal. An antenna, 
modeled as a filter, will have a similar effect. Some UWB 
antennas are fairly linear, such as Vivaldi antennas. However, 
other antennas, such as the LPA, have non-linear group delays 
(due to the frequency dependent phase-center of the LPA). 
A non-linear group delay is indicative of an ability for the 
antenna to store energy and causes resonances [12]. These 
resonances will result in oscillations and “ringing” of the 
output time signal. When used for a communications link or 
similar purposes, the resonance is a hindrance to performance 
and antenna engineers frequently try to mitigate or linearize 
the antenna group delay. However, for classification purposes, 
the frequency dependent resonance adds to the character of 
the antenna and can result in improved classification. The 
group delay is based on the phase difference between the 
input and output signals, which is directly affected by the 
antenna loading conditions, therefore the group delay will 
assist in classifying different termination conditions. Addition- 
ally, the bandwidth of the interrogation waveform and antenna 
frequency response will have a direct effect on group delay. 
A narrow-band interrogation pulse will experience a constant 
linear group delay whereas a wider bandwidth will experience 
the non-linear group delays, which will translate to a distorted 
signal that is useful for classification purposes. 


IV. ANALYTICAL MODEL DERIVATION 


The relationship between the received voltage of an antenna, 
Urx(t, 7,0, Y) and the electrical field incident on the antenna 
boresight, e;(t,r,0,w), can be modeled in both the time 
domain and frequency domain as [14]: 


tettre, V) neoa alent Y) 6) 
ve V0 

Urson 90) _ ye. g y Biles t 9 Y) 
SE THe we (6) 


where Ze and Zo are the antenna characteristic impedance 
and free space impedance respectively. Generally, h(t, 0, Y) 
and H(w,0, Y) are dependent on the direction of arrival of 
the incident wave in (6,1) space, but since we are only 
concerned with the bore sight arrival, the model for h(t, 0, y) 
and H(w,0,~) can be limited for model simplicity. 
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Antenna transmission is similar to reception, but with a 
significant difference. The transmit impulse response is ac- 
tually the first-order time derivative of the receive impulse 
response [14]. This is a consequence of the well known 
Lorentz Reciprocity constraint applied to Maxwell’s curl equa- 
tions [15]. Therefore, the time domain and frequency domain 
transmission antenna response can be modeled as [14]: 
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where e; is the electric field intensity at distance r from 
the antenna and u;, is the voltage incident on the antenna 
input terminals. Eq. 7 has a convolution with the Dirac-delta 
function, which is necessary to account for slowing due to the 
speed of light, c. As indicated above, the impulse response for 
transmission is related to the time-derivative of the impulse 
response for the receive case. This means that the antenna is 
going to radiate the time-derivative of the input voltage, which 
in the proposed circuit model, is the reflected voltage off of 
the antenna load. The gain of the antenna, G(f,0, Y) can be 

determined by the following [12]: 
An 


GU,0, 4) = EË ay o,wyP ©) 
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Fig. 2 shows the overall transmission line model for both 
the transmit and receive modes of the antenna. The use of two 
models is necessary because of the time derivative function of 
the transmit signal inherent to UWB antennas. 

For the proposed transmission line model, the antenna trans- 
fer function, H (w, 6, ~), will be modeled as a band pass filter. 
Some texts, such as [16] use the notation A(k,w) and h(k,w) 
to denote the transfer function (referred to as the spatio- 
temporal transmit and spatio-temporal receive characteristic 
respectively) where È is the incident field propagation direction 
instead of defining it in (0, Y) space as indicated above. 

An ideal UWB antenna will have a smooth gain across the 
entire frequency band which will eliminate dispersion during 
transmission. However, this is not realizable. Because of the 
Lorenz reciprocity constraint, 2jwh(—k,w) = coA(k,w), 
which has a temporal derivative in the jw term, the term 
A(w) = 0 when w = 0, and therefore hinders antenna radiation 
near w = 0 [15]. Additionally, realization of a antenna with a 
relatively constant gain across the entire frequency spectrum 
is impossible, and as such all antennas are ultimately band 
limiting. 

A typical UWB antenna has a >6:1 frequency range (i.e., 
the higher frequency is six times greater than the lower 
frequency) [17]. This bandwidth results in significant issues 
pertaining to linearity and impedance matching across the 
operating frequency of the antenna [16]. For the purposes 
of the transmission line model and simulations, a perfect 
impedance match between the antenna and transmission line, 
prior to termination, is assumed. However, further anaylsis 





can be performed to refine the model to include impedance 
mismatches across the antenna operating frequency. 


A. Re-radiation 


The voltage induced in the received antenna by the incident 
electric field generates a current through both the antenna 
and load impedances of Fig. 2. Part of the signal power is 
transferred into the transmission line, and part of the power is 
absorbed and re-radiated by the antenna impedance. Secondly, 
part of the signal is reflected off of the termination impedance 
which is also re-radiated by the antenna. Both of these reverse 
traveling waves can constructively or destructively interfere 
with the forward propagating signal (V,*). Assuming that the 
antenna is perfectly matched to the transmission line (unlikely, 
but reasonable for modeling) than the magnitude of the first 
reverse traveling wave is minimized leaving only the reflection 
off of the termination, (V ). 

The power delivered by the antenna in receive mode at a 
location rı is (assuming matched polarizations) [15] 
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where b’ is the source strength impacting the antenna in 
receive mode defined as [15] 
bps ile 
(w) = VaTpin kw) Ee 
0 
The process of illuminating the antenna with an incident 
plane wave results in an antenna-mode scattering process that 
consists of reception-reflection-reradiation [18]. This can be 
modeled in the frequency domain using the antenna transfer 
functions previously identified. A significant complication 
though is that the reflection coefficient of the termination 
is not generally constant across the frequency range, as is 
typical of a narrow-band system [19]. The calibration “short” 
and calibration “open” loads that have been experimentally 
examined generally have a constant reflection coefficient 
across the UWB frequency range. However, when dealing with 
other loads, the reflection coefficient will have a frequency 
dependence. The complete derivation of the resultant E-field 
at a distance, r, from the target antenna is derived as follows. 
The channel transfer function can be modeled as Eq. 12 
[12] where r is the distance of the incident propagating inter- 
rogation wave from the target antenna as well as the distance 
the resultant re-radiated electromagnetic field is measured and 
used for the noise radar cross-correlation and co is the speed 
of light. 


a1) 


ejwr/co 
m= (12) 
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Using Eq. 6, the voltage at the antenna terminals, U gy, due 
to the incident TEM wave propagating in the modeled channel, 
H., (V* in Fig. 2), is derived as: 


Ursl(w) = Hre(w, Ore, Vre) He(w)- V Zerz/Zo =V+(w) 
(13) 
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Because the distance the forward propagating TEM wave 
travels along the transmission line axis, z, is equivalent to the 
distance the reflected reverse propagating wave will travel, the 
reflected TEM wave can be derived simply as: 


V-(w) = V+ (w) -T(w) (14) 


In which I'(w) is the frequency dependent reflection co- 
efficient of the load, Zz. The reflected TEM wave is the 
voltage incident on the antenna terminals for re-radiation, 
Ur, = Ur, -T(w). Thus, the re-radiated field strength, 
Ers(w,r) at a distance r from is: 


Er,(w,r) = He (w) Hra (w, Ore, Yr2)jwU re (w) V Zo/Zere 
(15) 
Assuming that the characteristic impedance of the an- 
tenna, Ze is the same in both transmit and receive, as such 
Zere = ZeTz = Ze. Also, the transfer function of the 
antenna is the reciprocal in both transmit and receive with 
the exception of a sign change indicative of reverse signal 
flow. Therefore, Hr, = Hrs = Ha, and the overall equation 
relating the incident propagating TEM wave to the re-radiated 
electromagnetic field is: 


Erz(w,r) = H.( Harz, Orr. Vre) Arr (w, Ore, Ure) 
A, (w) - jwE;(w,r, 0, Y) Tw) (16) 
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B. Full System Characterization 


A plane wave propagating towards the bore sight of the 
antenna, while absolutely necessary for understanding the 
physics of the target-antenna model, is not practical for use 
with the RF-DNA tools since a second antenna is necessary 
for reception of the transmitted signal reflected off the antenna 
load, Z;. The two antennas therefore form a reciprocal link 
that can be fully characterized in the frequency domain using 
a similar transfer function approach. Assuming that both 
the transmit and receive (target) antennas are identical and 
operating in their respective far fields at a distance r, and that 
the bore sights of each antenna are facing each other, then the 
total analytical description of the transmitting propagation link 
in the frequency domain is as follows. 
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Note that 18b is derived from the reflection coefficient 
of a frequency dependent boundary, V,(f) = Vi(f)P'(f) as 
indicated in [20]. Combining Eqs. 18, and assuming that the 


antennas are matched, results in the following full system 
derivation: 
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Again, assuming Zerr = ZeTz = Ze and Hrs = Hrs = 
H,, the model reduces to: 
2jwre 


=’ 20 
Garo Uraa (f) (20) 

The time domain transmit and receive voltages, uTs(t) 
and up,(t) can be determined using the Inverse Fast Fourier 
Transform (IFFT). 


2 € 
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V. SIMULATIONS 


The system response will be simulated using Matlab. In 
order to properly characterize the impact of bandwidth on 
the RF-DNA classification process, the full radar link chain 
will be simulated. Therefore, transfer functions for the DNR 
transmit and receive antennas, Hyp, will be developed with 
an approximate bandwidth of 400-750 MHz as described in 
[5], with the assumption that matched antennas are used. 
Secondly, the target antenna will be characterized by a separate 
transfer function, Hr. Simulating the antenna as a filter will 
produce the desired frequency response, but will not simulate 
the dispersive nature of the antenna. Instead of characterizing 
the system in the frequency domain using multiple transfer 
functions, the system will be characterized in the time domain 
using Matlab’s filter command and a series of cascaded 
digital filters. This allows for easier control the frequency 
response of the DNR transmit and receive antennas as well 
as the response of the target LPA antenna. A digital all-pass 
filter with an arbitrary group delay will be implemented to 
produce the desired time dispersion. 

The input noise sequence is produced using Matlab’s 
randn function which creates an array of simulated random 
White gaussian Noise (WGN). In order to restrict the band- 
width of the signal, a 10th order digital Butterworth filter using 
the same butter and filter commands will be used to 
create the equivalent filter model for the UWB antennas. Both 
the center frequency of the WGN noise signal and bandwidth 
can easily be adjusted to accommodate various illumination 
signals. 

A final assumption made by this model, is to not include 
the effects of structural mode antenna scattering. It is not 
uncommon for an antenna structure to scatter more power than 
it absorbs [21]. However, that’s assuming scatter in a sphere 
radiating from the antenna whereas the amount propagating 
directly back towards the interrogation signal source can ac- 
tually be quite small. Scattering is also bandwidth dependent, 
since the power delivered to the load impedance is a function 
of the antenna input impedance which can vary with frequency. 
However, the structural mode scatter is heavily dependent on 
antenna geometry and location relative to other objects (such 
as a ground plane). The time-domain response of the structural 
mode scatter requires a more thorough understanding of the 
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target antenna and is heavily rooted in numerical solutions 
to Maxwell’s equations, which would generally require the 
use of computer programs such as CST. While the structural 
mode scatter will add to the character of the antenna, its 
necessarily complicates the issue of how bandwidth will 
affect the RF-DNA fingerprints and associated classification 
algorithms. 


A. Model Correlation Performance 


The received signal from the full system model is correlated 
with Matlab’s xcorr function. The results of the xcorr 
function applied to return signals off a modeled LPA with 
various noise interrogation signal bandiwidths is shown in 
Fig. 4. These results for the 100 MHz and 1 Ghz signal 
observationally correlate well with the correlation responses 
during previous experimentation [1], [3]. 


B. Antenna Type Classification 


For the first simulation, four antennas are iteratively interro- 
gated. The center frequencies of the antenna-filter bandwidths 
are spaced at 450, 525, 600 and 675 MHz with a band- 
width of 750KHz each. The interrogation signal is centered 
at fe = 562M Hz, or approximately in between the 2nd 
and 3rd antenna relative frequency spacing, and the signal 
bandwidth is increased from 100KHz to 1GHz. Fig. 5 shows 
the results of this interrogation, as a percent of the MDA/ML 
classifier tool properly classifying the validation fingerprint 
against the “true” class that fingerprint belonged to. A new 
classification model is constructed at each step increase in 
bandwidth. This approach was considered optimal because the 
user controls the interrogation waveform, and can choose the 
bandwidth that results in the optimal model performance. As is 
evident, none of the bandwidths achieved a classification above 
50%. This is because of the relative location of the center 
frequency of the interrogation signal. As the interrogation 
signal bandwidth is widened, it excites the passband of the 2nd 
and 3rd antennas, but since they are identical antennas, they 
confuse with each other. As the bandwidth is increased, the 
Ist and 4th antennas are also excited, but again, they confuse 
with the response of the 2nd and 3rd antennas resulting in total 
confusion. Changing the center frequency of the interrogation 
signal to identically match the center frequency of a specific 
antenna significantly increased the classification performance 
of that antenna. This is shown in Fig. 6. In this simulation, 
the antennas remained the same, but the center frequency 
of the interrogation signal was increased to 600MHz. The 
3rd antenna, with the bandwidth of 600-600.75MHz classified 
at above 90% across the span of interrogation bandwidths, 
whereas the other 3 antennas all remained at less than 40%. 
For this simulation, the bandwidth of the transmit antennas 
was set to approximately mirror the bandwidth of the DNR of 
375-750MHz. 

The second simulation extended the previous results, but 
with wideband antennas. Again, the antennas are spaced 
evenly with the previous center frequencies, but the bandwidth 
is increased to 150MHz for each antenna. The results are 

















Fig. 5. MDA/ML classifier performance for 4 broadband antennas with an 
interrogation signal center frequency of 562MHz 











Fig. 6. MDA/ML classifier performance for 4 broadband antennas with an 
interrogation signal center frequency of 600MHz 


shown in Fig. 7. The 3rd antenna (centered at 600MHz - red 
in the chart) classified at over 90% across the interrogation 
signal iterative bandwidth, as expected based on the previous 
narrower band simulation. The 2nd and 4th antennas (450- 
600MHz and 600-750MHz bandwidth respectively) did not 
classify as well, and were often confused with each other as 
opposed to the 3rd antenna. The first antenna was surprising, 
achieving 100% across the board. Analysis of the resultant 
confusion matrices indicated that it did not confuse with any of 
the other three antennas. The filter model used to represent the 
antenna was a 4th order Butterworth filter implemented using 
Matlab’s butter and filter commands. The simulated 
filter response only has about a 10dB magnitude difference 
between the peak of the passband and the center frequency 
of the interrogation signal (600MHz). xcorr is still able to 
derive a correlation result that is starkly different from the 
other three due to the center frequency spacing and band- 
width. Rerunning the simulation with an increased filter order 
(N=10), results in significantly improved classification for all 
4 antennas, as shown in Fig. 8. Additionally, the classification 
performance is significantly less responsive to the interrogation 
signal bandwidth. Based on this, more important than the 
interrogation signal bandwidth is the frequency response of 
the target antennas. 


C. Signal Attenuation Classification 


In this simulation, the ability of the RF-DNA fingerprints 
with classifier algorithms to classify an attenuator placed 
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Fig. 3. Measured versus Simulated Impulse Response of the Transmit and Receive DNR LPA antennas with simulated Group Delay 











Fig. 7. MDA/ML classifier performance for 4 wideband antennas with an 
interrogation signal center frequency of 600MHz (4th order filter approxima- 
tion) 














Fig. 8. MDA/ML classifier performance for 4 wideband antennas with an 
interrogation signal center frequency of 600MHz (10th order filter approxi- 
mation) 


between the antenna and the load is qualitative analyzed. To 
simulate the effect of the attenuator, the signal both before 
and after reflection from the load is attenuated by a scaling 
factor. Five different scale factors are used equivalent to the 





following values for the attenuator, ag = [0 3 6 9 12] dB. 
Initially, the DNR antennas and target antennas are assumed 
identical with a bandwidth of 375-750MHz. The noise signal 
has a center frequency, fe = 562 MHz, in the middle of the 
antenna bandwidths. The bandwidth of the noise is varied 
from a bandwidth of approximately 90KHz to 900MHz. It is 
anticipated that the increased bandwidth will improve the per- 
formance of the classifier since the peak value of the sinc will 
scale larger with increased bandwidth, resulting in an increased 
instantaneous amplitude response of the RF-DNA fingerprints 
in the correlation region. Also, since the noise signal is filtered 
to produce the desired bandwidth of the interrogation signal, 
the amount of energy contained within the instantaneous 
bandwidth is increasing as the bandwidth increases. Fig. 9 
shows the percentage of time that the MDA/ML classifier 
correctly classified each attenuator value during model val- 
idation. Again, a new classification model is constructed for 
each step increase in interrogation signal bandwidth. The chart 
indicates that a narrow band signal has difficulty classifying all 
5 attenuator values. As the bandwidth is increased, the ability 
to identify the first two attenuator values increased quickly, 
and the third increases as well, although not as fast. The 
9dB and 12dB attenuation values do not increase significantly, 
probably because the correlation response is too small for 
those attenuator values and they confuse with each other 
across the interrogation signal bandwidth. These simulation 
results indicate interrogation signal bandwidth directly affects 
RF-DNA classification performance. 


VI. CONCLUSIONS 


In conclusion, a literature review was conducted on the use 
of filters as antenna models, and the appropriate characteristics 
(such as the impulse repose, group delay and scatter) were 
evaluated. An analytical expression for the transfer func- 
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Fig. 4. Simulated XCORR System response to a <1MHz (top left), 1OMHz (top right), 100MHz (bottom left) and 1GHz (bottom right) bandwidth interrogation 
signal 


filter model. Overall, the interrogation signal instantaneous 
bandwidth has a significant impact on classifier performance. 
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